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Abstract 

We find using Monte Carlo simulation the phase structure of noncommutative U(l) 
gauge theory in two dimensions with the fuzzy sphere S% as a non-perturbative regulator. 
There are three phases of the model, i) A matrix phase where the theory is essentially 
SU(N) Yang-Mills reduced to zero dimension . ii) A weak coupling fuzzy sphere phase 
with constant specific heat and Hi) A strong coupling fuzzy sphere phase with non-constant 
specific heat. The order parameter distinguishing the matrix phase from the sphere phase 
is the radius of the fuzzy sphere. The three phases meet at a triple point. We also give 
the theoretical one-loop and i expansion predictions for the transition lines which are in 
good agreement with the numerical data. A Monte Carlo measurement of the triple point 
is also given. 
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1 Introduction 

Quantum noncommutative ( NC ) gauge theory is essentially unknown beyond one-loop 
PP. In the one-loop approximation of the quantum theory we know for example that gauge 
models on the Moyal- Weyl spaces are renormalizable These models were also shown to 
behave in a variety of novel ways as compared with their commutative counterparts. There 
are potential problems with unitarity and causality when time is noncommuting, and most 
notably we mention the notorious UV-IR mixing phenomena which is a generic property of all 
quantum field theories on Moyal- Weyl spaces and on noncommutative spaces in general |T]|3]. 
However a non-perturbative study of pure two dimensional noncommutative gauge theory was 
then performed in 0. For scalar field theory on the Moyal- Weyl space some interesting non- 
perturbative results using theoretical and Monte Carlo methods were obtained for example in 
jSj. An extensive list of references on these issues can be found in |T] and also in 0] 

The fuzzy sphere ( and any fuzzy space in general ) provides a regularized field theory in the 
non-perturbative regime ideal for Monte-Carlo simulations. This is the point of view advocated 
in [7j. See also [E1III1EI] for quantum gravity, string theory or other different motivations. These 
fuzzy spaces consist in replacing continuos manifolds by matrix algebras and as a consequence 
the resulting field theory will only have a finite number of degrees of freedom. The claim is that 
this method has the advantage -in contrast with lattice- of preserving all continous symmetries 
of the original action at least at the classical level. This proposal was applied to the scalar 
4 model in jTTj. Quantum field theory on fuzzy spaces was also studied perturbatively quite 
extensively. See for example [T3*l ITU IT5] . For some other non-perturbative ( theoretical or 
Monte Carlo ) treatement of these field theories see . 
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Another motivation for using the fuzzy sphere is the following. The Moyal-Weyl NC space 
is an infinite dimensional matrix model and not a continuum manifold and as a consequence it 
should be regularized by a finite dimensional matrix model. In 2 dimensions the most natural 
candidate is the fuzzy sphere S% which is a finite dimensional matrix model which reduces 
to the NC plane in some appropriate large iV flattening limit. This limit was investigated 
perturbatively in P^ IT7j for scalar and Yang- Mills field theories respectively. In 4— dimensions 
we should instead consider Cartesian products of the fuzzy sphere [TH], fuzzy CP^ [TH] or 
fuzzy S 4 An alternative way of regularizing gauge theories on the Moyal-Weyl NC space 
is based on the matrix model formulation of the twisted Eguchi-Kawai model. See for example 

The goal of this article and others ^2JI22J is to find the phase structure ( i.e map the different 
regions of the phase diagram ) of noncommutative U(l) gauge theories in 2 dimensions on the 
fuzzy sphere S%. There are reasons to believe that the phase diagram of NC U(N) models 
will be the same as that of their U(l) counterparts thus we will only concentrate on the U(l) 
models. Furthermore it seems that the nature of the underlying NC space is irrelevant. In other 
words U(l) gauge models on the NC Moyal-Weyl plane Rj , on the fuzzy sphere and on 
the NC torus T# will fall into the same universality class. Hence we consider solely the fuzzy 
sphere since it is the most convenient two dimensional space for numerical simulation. 

There seems to exist three different phases of U(l) gauge theory on S^. In the matrix phase 
the fuzzy sphere vacuum collapses under quantum fluctuations and there is no underlying sphere 
in the continuum large N limit. Rather we have a U(N) YM theory dimensionally reduced to 
a point. In this phase the model should be described by a pure ( possibly a one-)matrix model 
without any spacetime or gauge theory interpretation. This phenomena was first observed in 
Monte Carlo simulation in j2H| for m — 0. In [22] it was shown that the fuzzy sphere vacuum 
becomes more stable as the mass m of the scalar normal component of the gauge field increases. 
Hence this vacuum becomes completely stable when this normal scalar field is projected out 
from the model. This is what we observe in our Monte Carlo simulation in the limit m — >oc. 

The principal new discovery of this paper is that the fuzzy sphere phase splits into two 
distinct regions corresponding to the weak and strong coupling phases of the gauge field. These 
are separated by a third order phase transition. This transition is consistent with that of a 
one-plaquette model j2H|. Our results indicate that non-perturbative effects play a significant 
role than expected from the 1/N study of [24 . In particular these results indicate that quantum 
noncommutative gauge theory is essentially equivalent to ( some ) quantum commutative gauge 
theory not necessarily of the same rank. This prediction goes also in line with the powerful 
classical concept of Morita equivalence between NC and commutative gauge theories on the 

torus mini- 

This article is organized as follows. In section 2 we will describe the phase diagram of the NC 
£7(1) gauge model in 2D. In section 3 we will review the one-loop theory of the model. In section 
4 we will discuss our Monte Carlo results with some more detail. In section 5 we will introduce 
the one-plaquette approximation of the model and then we will give a theoretical derivation of 
the one-plaquette line. We conclude in section 6 with a summary and some general remarks. 
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In the appendix we discuss ( among other things ) the measurement of order parameters and 
probability distribution. 



2 Phase diagram 

The basic action is written in terms of three NxN matrices X a as follows 



1 „ r „ 2ia 



-Tr[X a ,X b ] 2 + ^e abc TrX a X b X c 



Nm 
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Nm 2 a 2 TrX 2 + Tr(X 2 y. (1) 

/Co 



S = N 

The basic parameters of the model are a = ay/N and m. The gauge coupling constant is 



g 2 = -Jj. We will also need a = aJl — j|, a = ay/N and rh = — . The one-loop critical 
values of 5 , a and a are 5*, d* and a* respectively. Clearly in the large N limit d* = d* 
and a* = a^y/N . For reasons which will become clear in the text the measured critical values 
are denoted as follows. The measurement of a* is denoted by a s . There are two physically 
distinct measurements of a* denoted by a ma and a mi . In terms of a these are given by a ma = 
a ma \J 1 — jj and a m i = a m iJl — j^. There are two physically different measurements of a* 
denoted by ttp — a ma y/ N and a s — a s y/ N . 

The phase diagram of the model (JIJ is given in figure (1). This is the central result of this 
article. In this section we will briefly explain the main properties of the different phases of the 
model. More detail will be given in the rest of the article. 

We measure the average value of the action < S > as a function of a and we measure the 
specific heat C v =< S 2 > — < S > 2 as a function of a for different values of N. We consider 
N = 4,6, 8, 10, 12, 16. In the first step of the simulation the mass parameter m is taken to 
be some fixed number. Then we vary the mass parameter and repeat the same experiment. 
The choice of a for the specific heat is only due to finite size effects and has no other physical 
significance since in the large N limit a — a. 

We observe that different actions < S > which correspond to different values iV ( for some 
fixed value of m ) intersect at some value of the coupling constant a which we denote a s . In 
the limit of small masses, viz m — >0, this intersection point marks a discontinuity in the action 
and it occurs around the value a s = 2.2. In figure 2 we plot the action < S > versus a for 
N = 10, 12, 16 and m 2 = 0.25, 3, 100. For large masses we observe that the intersection point 
becomes smoother. It is also clear that the critical point a s decreases as we increase m and it 
will reach in the limit m — >oo. 

For the specific heat the situation is more involved. We observe in the limit m — >0 a peak 
around a s = 2.2 which marks a sharp discontinuity in C v . See the first graph in figure 3 or 
figure 7. Above this critical value the specific heat is given by C v = N 2 while below this critical 
value it is given by C v = 0.75N 2 . The regime a>a s is the fuzzy sphere phase whereas a<a s 
corresponds to the so-called matrix phase. 

As m increases things get more complicated and they only simplify again when we reach 
large values of m. a ma and a m i are precisely the values of a at the maximum and minimum 
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The phase diagram 
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Figure 1: The phase diagram of the model (JTJ). The two fits © and (0) are expected to coincide 
very well with the data only for very large masses. The fits (0) and (@J) are identical for large 
masses ( or equivalently past the triple point ). Above the upper critical line we have a fuzzy 
sphere in the weak regime of the gauge theory. Between the two lines we have a fuzzy sphere in 
the strong regime of the gauge theory. Thus the upper critical line is the one-plaquette critical 
line ©. Below the lower critical line we have the matrix phase. This last line agrees very well 
with the one- loop prediction (jSJ). The one-plaquette line approaches in the limit m — >oo a 
constant value given by loga p = log(3.35) = 1.21. The triple point is also seen to exist within 
the estimated range. Before we reach the triple point the critical line agrees as well with the 
one-loop prediction (JU) . Recall that log m 2 = log m 2 — 2 log N, log a p = log a ma + 0.5 log iV and 
log a s + 0.5 log iV = log a s . 
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values of the specific heat. Thus a ma and a m i are the values of a at the maximum and minimum 
values of the specific heat. The peak of the specific heat moves slowly to smaller values of at 
as we increase m. The agreement between a ma and a s for small masses is good whereas the 
values a m i at the minimum of C v for small masses are significantly different from a s . Thus in 
this regime of small masses a ma is still detecting the S^— to- matrix phase transition. Similarly 
to the case m = the specific heat C v as a function of a is equal to N 2 in the fuzzy sphere 
phase for values of a such that a>a ma . 

The physics is drastically different for large masses since the roles of a ma and a m i are 
completely reversed. There is a shallow valley in the specific heat starting to appear for values 
of a inside the matrix phase as m slowly increases. Furthermore as m keeps increasing we 
observe that the peak flattens slowly and disappears altogether when the mass becomes of the 
order of m 2 ~10. At this stage the well in the specific heat becomes on the other hand deeper 
and more pronounced and its minimum a mi is moving slowly to smaller values of the coupling 
constant a. By inspection of the data we can see that a m i and a s starts to agree for larger 
masses and thus a m i captures exactly the S^-to-matrix phase transition in this regime. The 
physical meaning of the the critical point a ma becomes also different for large masses where it 
becomes significantly different from a s . Since there is no peak the definition of a ma becomes 
different 1 . a ma is now the value of a at which the specific heat jumps and becomes equal 
to iV 2 . This is where the one-plaquette phase transition between weak and strong regimes of 
gauge theory on the fuzzy sphere occurs. In figure 3 we plot the specific heat versus a for 
N = 10, 12, 16 and m 2 = 0.25, 3, 100. In particular remark how the shape of the specific heat 
changes with m. 

As it turns out we can predict the S^-to-matrix phase transition from the one-loop theory of 
the model To this end we consider the following background matrices D a = a<j)L a where 4> 
is the radius of the sphere and L a are the generators of SU(2) in the irreducible representation 
— ip^. Then we compute the one-loop effective potential in the background field method [22] or 
by using an RG method [T2] . One finds the result 



2^,4 r 



Vi- 



N a 



loop 



1 + m 



4 



m 
~2 



+ N 2 log + O(iV). 



(2) 



It is not difficult to check that the corresponding equation of motion of the potential (J2J) admits 
two real solutions where we can identify the one with the least energy with the actual radius 
of the sphere. This however is only true up to a certain value <5* of the coupling constant at 
where the quartic equation ceases to have any real solution and as a consequence the fuzzy 
sphere solution D a = a<pL a ceases to exist. In other words the potential below the value a* of 
the coupling constant becomes unbounded and the fuzzy sphere collapses. The critical values 
can be easily computed and one finds in the limit m — >0 the values 0* = 0.75 and a* = 2.09. 

1 It is not difficult to see that this new definition is consistent with the previous definition of a ma - 
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Extrapolating to large masses we obtain the scaling behaviour 



and 



7=2 (3) 



III* + V2-1 



In other words the phase transition happens each time at a smaller value of the coupling 
constant a and thus the fuzzy sphere is more stable. This one-loop result is compared to the 
non-perturbative results a s and a m i coming from the Monte Carlo simulation of the model (Q) 
in figure 4. As one can immediately see there is an excellent agreement between the three values 
in the regime of large masses as discussed above, a* and a s agree as well for small masses. See 
also the phase diagram (1). 

For large values of m the scaling of the coupling constant a as well as of the mass parameter 
m is found to differ considerably from the m = case. It is now given by a = ay/~N, m = ^ 
The above theoretical fit (jlj will read in terms of a and fh as follows 

= ( 5 ) 
m* 

This is the fit used for the lower critical line in the phase diagram (1). The critical value 
a s = a s yN falls nicely on the top of this fit for all values of the mass. 

The fit of the critical value a ma for m small is given by equation (j3J). Thus we expect 
agreement between a* and a s from one hand and a ma from the other hand in the range of 
small masses. For m large we find that we can fit the data a ma to ( recall that a p = a ma \fN ) 

rO.04.1 

a p = 3.35±0.25+ — - R 6 

m 2 

In other words in the limit m — >oo we can fit the data to the line a p = 3.35. This is what we 
call the one-plaquette critical line. See figure 5. This is the fit used for the upper critical line 
in the phase diagram (1). 



3 The one-loop calculation 



In this section we will follow 

We are interested in the most general gauge theory up to quartic power in the gauge field 



on the fuzzy sphere S 2 L+l . This is obtained as follows. Let X a 



1,2,3, be three iVxiV 



hermitian matrices and let us consider the action 



-Tr[X a ,X b f + 2 -^-e abc TrX a X b X L 



2\2 



+ (3TrX^ + MTr(X>) 



(7) 
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The action for m =0.25 and N=10,12 
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The action for m 2 =100 and N=10,12,16 




jure 2: The action for m 2 = 0.25, 3, 100 and N = 10, 12, 
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The specific heat for m =0.25 and N=1 0,1 2 
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Figure 5: The one-plaquette critical line. 
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This action is invariant under the unitary transformations X a — >UX a U + . This model is also 
invariant under SU{2) rotational symmetries X a — >gX a g + = R ab (g)X b where the group ele- 
ment g is given explicitly by g = exp(iu a L a ) for some constant vector uo and R(g) is the spin one 
irreducible representation of g. L a are the generators of 577(2) in the irreducible representation 
2". They satisfy [L a ,L b ] = ie abc L c , L 2 = c 2 = \{k + 1) and they are of size (L + 1) x (L + 1). 
This action is bounded from below for all positive values of M and the trace is normalized such 
that Trl = N = L + 1. 

The a, j3 and M are the parameters of the model. We are interested in the particular case 
where f3 = —2Mc 2 a 2 . In this case the potential becomes 

V = (3TrXl + MTr(Xl) 2 = MTr{X 2 a - c 2 a 2 ) 2 - MNa 4 c 2 2 . (8) 

In most of this article we will discuss U(l) gauge theory on the fuzzy sphere. We start with the 
values M — (5 — 0. This corresponds to the Alekseev, Recknagel, Schomerus action obtained in 
effective string theory describing the dynamics of open strings moving in a curved background 
with S 3 metric in the presence of a Neveu-Schwarz B-field. We notice that with M = (3 = the 
trace part of X a simply decouples. As a consequence we can take X a to be traceless without 
any loss of generality. The classical absolute minimum of the model is given by 

X a = aL a (9) 

where L a are the generators of 577(2) in the irreducible representation ^=^-^-.The quantum 
minimum is found by considering the configuration X a = a<j)L a where the order parameter a<ft 
plays the role of the radius of the sphere with a classical value equal a. The complete one-loop 
effective potential in this configuration is given in the large N limit by the formula ( with 

a = VNa ) 



^■1 



V&(4>) = 2c 2 a 



1 jA ^3 



+ 4c2 log0 + subleading terms. (10) 



It is not difficult to check that the equations of motion admits two real solutions where we can 
identify the one with the least energy with the actual radius of the sphere. However this is only 
true up to a certain value a* of the coupling constant a where the quartic equation ceases to 
have any real solution and as a consequence the fuzzy sphere solution @ ceases to exist. In 
other words the potential V e s below the value a* of the coupling constant becomes unbounded 
and the fuzzy sphere collapses. The critical value can be easily computed and one finds 

0* = - , 5* = 2.08677944. (11) 

Now we add the potential term (jSJ) with mass parameter 2M = Nm 2 /c 2 . In this case the 
matrices X a can not be taken traceless. The effective potential becomes 



V eff = 2c 2 a 



~0 4 -~0 3 + ~™ 2 (0 2 -l) 2 
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+ 4c 2 log + ^Tr 3 TRlogA. (12) 



TR is the trace over 4 indices corresponding to the left and right actions of operators on matrices 
of size L + l while Tr 3 is the trace associated with the action of 3— dimensional rotations. The 
Laplacian A in the gauge = 1 + — is given by 



A 



C 2 + (- - 1)(J 2 - C 2 - 2) + 2m 2 (l - — ). 



(13) 



The eigenvalues of C 2 ( which is the Laplacian on the sphere ) and J 2 ( which is the total angular 
momentum on the sphere ) are given respectively by 1(1 + 1) and j(j + 1) where I = 1, L and 
j —1 + 1,1,1 — 1. The corresponding eigentensors are the vector spherical harmonics operators. 
Let us also notice that from the requirement that the spectrum of A must be positive we can 
derive a lower and upper bounds on the possible values which the field <ft can take. For example 
for m 2 = we can find that 2/3 < <fi < 3. 

We can show ( at least ) for small values of the mass m that the logarithm of A is subleading 
in the large N limit compared to the other terms and thus the potential reduces to the simpler 
form 

I .„ 1 1 



V, 



cff 



2c 2 a 4 



4 v 



1 



+ 4c 2 log . 



(14) 



Solving for the critical value using the same method outlined previously yields the results 



8(1 +m 2 ) 



1 + 




32m 2 (1 ■ 
9~ 



m 



2\i 



(15) 



(16) 



Extrapolating to large masses (m — >oo) we obtain the scaling behaviour 



a* 



m- 



+ V2-1 



(17) 



In other words the phase transition happens each time at a smaller value of the coupling 
constant a and thus the fuzzy sphere is more stable. 

It is therefore sensible to expand the action (jjj) around the fuzzy sphere solution (0) by 
introducing a U(l) gauge field A a on the fuzzy sphere S% as follows X a = a(L a + A a ). The 
action becomes 



S 



N 



Ag 2 N 
1 



Tr 



F^ + i[A a ,A b 



1 



2# 2 iV 



e„hrTr 



l -F { ^A c + l -[A a ,A b }A c 



+ 



2m 2 

7n 



7W 



~4 ± ~4 2 

-a c 2 - -at c 2 m . 
o 2 



$ is the normal covariant scalar component of the gauge field on the fuzzy sphere defined by 
y/4c2& = L a A a + A a L a + A 2 . F a b = + zL4 a ,Ab] is the U(l) covariant curvature where 
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Fab = i[L a , Ab] — i[Lb, A a ] + e a bcA c and g is the gauge coupling constant denned by \ = a 4 . In 
the continuum limit L — >oo all commutators vanish and we get a U(l) gauge field coupled to 
a scalar mode $ = ft. A with curvature = iC a A b — iCbA a + e a b c A c where C a = —i^abc n b^ ■ 
We find ( by also neglecting the constant term ) 

s ~ - H^-^-J^^i^- (19) 

The quantization of the fuzzy action Sn yields a non-trivial effective action in the continuum 
limit which for generic values of the mass parameter m is different from S^. For example we 
have established by explicit calculation of the quadratic effective action the existence of a gauge- 
invariant UV-IR mixing problem in U (1) gauge theory on the fuzzy sphere for the value m — 0. 
We find 

^ - £ /, S^ 1 + 2 ^ - sW, + 2c/A3)A + '^L 

+ other non local quadratic terms. (20) 

The operator A 3 is a function of the Laplacian C? with eigenvalues A 3 (/c) given by k(k + 
l)A 3 (fc) = Ep =2 I- Clearly Tg^S^ which is the signature of the UV-IR mixing in this model. It 
is expected that the same result will also hold for generic values of the mass parameter m. 

The calculation can also be done quite easily in the limit m — >oo and one finds that there 
is no UV-IR mixing in the model in this case. The UV-IR mixing is thus confined to the 
scalar sector of the model since the limit m — >oo projects out the scalar fluctuation $. It 
is hence natural to think that the extra matrix phase observed in the phase structure of the 
theory is related to this mixing; in other words it is the non-perturbative manifestation of the 
perturbative UV-IR mixing property since as we have shown this phase seems also to disappear 
in the limit of large masses. 

4 Monte Carlo simulation 

From the above one-loop argument it is expected to observe at least one phase transition 
on the line m — f3 — 0. This is a continuous first order phase transition from a fuzzy sphere 
phase with a > a* to a pure matrix phase with a < a*. It is also expected that this critical 
value a* decreases with the value of m ( keeping j3 fixed equal — 2Mc2<y 2 with 2M = Nm 2 /c2 
) and it becomes zero when we let m — >oo. 

4.1 Zero mass 

We start with M — (3 — 0. To detect the different phases of the model we propose to measure 
the following observables. First we measure the average value of the action, viz < S >. Second 
the specific heat will allow us to demarcate the boundary between the different phases. It is 
defined by C v =< S 2 > - < S > 2 . 
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In order to determine the critical point (if any) we run several simulations with different 
values of N, say N = 4,6, 8, 10, 12, 16. We always start from a random (hot) initial configuration 
and run the metropolis algorithm for Ttherm + TcorrxTmont Monte Carlo steps. Ttherm is 
thermalization time while Tmont is the actual number of Monte Carlo steps. Two consecutive 
Monte Carlo times are separated by Tcorr sweeps to reduce auto-correlation time. In every 
step ( sweep ) we go through each entry of the three matrices Xi, X 2 and X 3 and update it 
according to the Boltzman weight. This by definition is one unit of time ( Monte Carlo time ) 
in the generated dynamics. For every N and a we tune appropriately Ttherm, Tmont, Tcorr 
as well as the interval / from which we choose the variation of the entries of the matrices X a 
so that to reduce auto-correlation times and statistical errors. 

The continuum limit of a given observable will be obtained by collapsing the corresponding 
data, in other words finding the scaling of this operator in the large N limit which yields an 
iV— independent quantity. For example the scaling of the coupling constant a with N is clearly 
given by a = \/N a as anticipated from the one-loop calculation. 

For the action the data is plotted in figure 6. We remark that the 4 curves with N = 4, 6, 8 
and 10 all intersect around the point 

a s = 2.2±0.1. (21) 

This is the critical point since it is independent of N as it should be. The collapse of the data 
is given by < S > j4ci- Indeed a very good fit for the action < S > is given by the classical 
action in the configuration X a = aL a , i.e 

< S >= (22) 

The data for the specific heat is shown on figure 7. We can immediately remark that C v peakes 
around the above critical point. More precisely the peak is at the values a = 2.25±0.05, 2.1±0.1 
and 2.1±0.1, 2.2±0.1 for N = 4,6 and 8, 10 respectively. 

From figure 7 the correct scaling of the specific heat is given by C v /4c2- Let us also remark 
that the specific heat is equal C v = N 2 — 1 in the fuzzy sphere phase and C v = 0.75(N 2 — 1) in 
the matrix phase. 

We have therefore established the existence of a first order phase transition from the fuzzy 
sphere to a matrix phase in agreement with the one-loop calculation. The next step is to add 
to the Alekseev-Recknagel-Schomerus model the potential term (JHJ). 

4.2 Non-zero mass : the S^— to-matrix phase transition 

As we have shown m? plays precisely the role of the mass parameter of the normal scalar 
field in the fuzzy sphere phase. From the one- loop calculation as well as from the large 1/N 
expansion it is argued that the fuzzy sphere becomes more dominant ( i.e it becomes more 
stable under quantum fluctuations ) as we increase the mass m of the scalar mode. In the 
limit m — >oo we expect the matrix phase to disappear altogether. In this limit m — >oo the 
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Figure 6: The action for zero mass for N 
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Figure 7: The specific heat for zero mass for N = 4, 6, 8, 10. 



normal scalar field decouples from the pure two dimensional gauge sector and as a consequence 
it is natural to conjecture that the matrix phase ( and correspondingly the perturbative UV-IR 
mixing phenomena ) has its origin in the coupling of this extra normal mode to the rest of the 
dynamics. Another way of putting this conjecture is that the presence of the matrix phase ( 
which is absent in the continuum theory ) is nothing else but a non-perturbative manifestation 
of the perturbative UV-IR mixing. However although this is true to a large extent there are 
more non-trivial things happening in this limit as we will now report. 

We again measure the action < S > for non-zero values of the mass m for N = 4, 6 and 8. 
The results are shown on figure 8. As before the action < S > is scaled as < S > /4c 2 as a 
function of a. It is also immediately clear that the critical point decreases as we increase m. In 
other words the fuzzy sphere becomes more dominant as promised by the one-loop calculation. 
For example for m 2 = 0.5 the actions for N = 4, N = 6 and N = 8 intersect at a s = 1.9±0.1. 
The theory predicts a critical value given by <5* = 1.72 which is reasonably close. 

We repeat the above calculation for various values of the mass m. The intersection point of 
the actions with different N defines the critical value a s . This value follows to a good accuracy 
the one-loop prediction given by equation (JTTj). As it turns out this phase transition is also 
captured by the minimum of the specific heat ( more on this below ). The phase diagram of 
the fuzzy sphere-to- matrix phase transition is shown on figure 9. A very good fit of < S > is 
given by the classical action in the fuzzy sphere configuration. For non-zero mass this is given 
by the expression 



< S > 



~4 

crc 2 
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m 
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The action for m =0.5 and N=4,6 and 8 
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$ure 8: The action for non-zero mass. The fit is given by equation 
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Figure 9: The phase diagram of the S 2 N — to- matrix phase transition. The fuzzy sphere phase is 
above the solid line while the matrix phase is below it. The critical point a s is the intersection 
point of the actions with N = 4, 6 and 8. This agrees with the one-loop prediction of the critical 
line ( the solid line ) given by equation (|17j) . For large masses a s coincides with the minimum 
a m i of C v . For small masses a s coincides with the maximum a ma of C v . 
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4.3 Specific heat: the one-plaquette phase transition 

As soon as the mass m takes a non-zero value the specific heat C v starts to behave in a very 
different way compared to its behaviour for zero mass. We observe a new phase transition for 
large enough masses which resembles very much the one-plaquette phase transition in ordinary 
2— dimensional gauge theory. This is measured by the maximum a ma of the specific heat. For 
small values of the mass parameter m the maximum a ma is defined by the position of the peak 
of C v . For large values of m the peak in C v disappears and a ma is given by the value of the 
coupling constant at which the specific heat discontinuously jumps to one. We will also measure 
the minimum a mi of C v which will capture the S^— to- matrix phase transition. 

For small values of m the scaling of the coupling constant a is found to differ only slightly 
from the m = case. It is given by 



In the large N limit a/ a — >1 and thus this different scaling is only due to finite size effects 
and has no other physical significance. This is expected for small masses. a ma and a m i are 
actually the values of a at the maximum and minimum of the specific heat. The peak of the 
specific heat moves slowly to smaller values of the coupling constant as we increase m. The 
agreement with the one-loop prediction given by equation (fTTj) as well as with a s is fairly good 
and thus in this regime a ma is still detecting the S^— to- matrix phase transition. Similarly to 
the case m = the specific heat C v /Ac2 as a function of a is equal to 1 in the fuzzy sphere 
phase. However there is a shallow valley starting to appear for values of a inside the matrix 
phase. The values a m i of the minimum of C v for these small masses are significantly different 
from a s ( see the phase diagram on figure 9 ). As an example the data for m 2 = 0.25,4.75 for 
N = 6, 8 is shown on figure 10. 

As m keeps increasing deviation from the one-loop prediction becomes important. The data 
for m 2 = 40, 200 for iV = 6, 8 is shown on figure 10. We observe that the peak flattens slowly 
and disappears altogether when the mass becomes of the order of m 2 ~10. 

Although the peak in C v disappears we know that the S^r-to-matrix phase transition is still 
present as indicated by the non-vanishing of a s ( from the phase diagram on figure 9 ) . The 
physical meaning of the the critical point a ma becomes different for large masses. This is where 
the one-plaquette phase transition between weak and strong regimes of gauge theory on the 
fuzzy sphere occurs. The valley in the specific heat becomes on the other hand deeper and 
more pronounced as m increases and its minimum a m i is moving slowly to smaller values of the 
coupling constant a. By inspection of the data ( phase diagram on figure 9 ) we can see that 
a m i and a s starts to agree for larger masses and thus a m i captures exactly the S^-to-matrix 
phase transition. 

The two regimes with m small and m large are thus physically distinct; in the first regime 
we have two phases : the fuzzy sphere phase a>a s and the matrix phase a<a s whereas in the 
second regime we have three phases. Beside the matrix phase for a<a m i we have two more 




(24) 
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phases where we have a stable fuzzy sphere as the underlying spacetime structure. These two 
phases correspond to U(l) gauge theory on the fuzzy sphere S% in the weak a>a ma and strong 
«mi<«<«ma regimes respectively. There exists therefore a triple point where the three phases 
coexist. 

For large values of m the scaling of the coupling constant a as well as of the mass parameter 
m is found to differ considerably from the m = case. It is now given by 

i — - m 

a = ay N , m = —. (25) 
The theoretical fit (jl7j) will read in terms of a and fh as follows 

a. = [4$. (26) 

We define the one-plaquette transition point by the value a ma of the coupling constant a ( or 
equivalently a in the large N limit ) at which the specific heat discontinuously jumps to one. 
In terms of a this is given at the value 

a p = a ma VN. (27) 

The fit of the critical value a ma for m small is given by equation (|17|) while for m large we find 
that we can fit the data to 

rO.04.1 

a v = 3.35±0.25+ — - 2 . 28 

In other words in the limit m — >oo we can fit the data to a p = 3.35±0.25. In the next section 
we will give a theoretical derivation of the value 3.35 from the one-plaquette approximation of 
gauge fields on the fuzzy sphere in the weak regime a>a p . 

These results are summarized in the phase diagram on figure 11. 

Finally we point out that we can estimate the values olt and fn\ of the coupling constant a 
and the mass parameter fh 2 at the triple point by equating the fits (J2T?|) and (|28|). We obtain the 
two solutions 1) mf. = 0.009 and ay = 5.46 or equivalently logm^ = —4.71 and loga^ = 1.7 
and 2) fh 2 -, = 0.001 and oit = 9.46 or equivalently logm^ = —6.91 and logar = 2.25. The 
triple point must therefore exist between these two points, viz 

1.7<loga T <2.25 (29) 

and 

-6.91<log7r4< -4.71. (30) 

The most important remark we can draw from this calculation is that the fuzzy sphere phase 
bifurcates into two distinct phases ( the weak coupling and the strong coupling phases of the 
gauge field ) almost as soon as we tune on a non-zero mass. The models with and without a 
mass term are indeed very different. 
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Figure 10: The specific heat for different values of m 2 
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The phase diagram ( logarithmic scale ) 
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Figure 11: The phase diagram of the one-plaquette phase transition. For large values of m 
the a p — line is the critical one-plaquette transition line while the a m i— line ( or equivalently 
the a s — line ) is the critical fuzzy sphere-to-matrix transition line. For small values of m the 
ot ma — line coincides with the a s — line which is the critical fuzzy sphere-to-matrix transition line 
in this regime. These two lines seem to bifurcate at the triple point. 
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5 The one-plaquette model and 1/N expansion 



In this section we will follow [T2*j . 



5.1 The one-plaquette variable W 

We start by making the observation that in the large m — >oo limit we can set $ = since 
the normal scalar field becomes infinitely heavy ( m is precisely its mass ) and thus decouples 
from the rest of the dynamics. Hence we can effectively impose the extra constraint X 2 = o?Ci 
on the field X a in this limit m — >oo. The action (J7J) with j3 = —2Mc<iO? becomes in the limit 
m — >oo first and then iV — >oo a commutative U(l) action on the ordinary sphere. 

The aim is to relate the action (JJJ) with the one-plaquette action. To this end we introduce 
the 2Nx2N idempotent 

7 = -^(l2tf + 2<7 L„) , 7 2 = 1 (31) 

where a a are the usual Pauli matrices. It has eigenvalues +1 and —1 with multiplicities N + 1 
and N — 1 respectively. We introduce the covariant derivative D a = L a + A a through a gauged 
idempotent 7^ as follows 



Id = 7" 



/ ^y2 



7 = 1(1 + 2a a D a ) = 7 + ^a a A a , f = 1 + + l_e abc a c F ab . (32) 

Since we are interested in the large m — >oo limit we set $ = 0. Clearly 7^ has the same 
spectrum as 7. In fact 7^ is an element of the d^— Grassmannian manifold U(2N)/U(N + 
l)xU(N — 1) and hence it contains the correct number of degrees of freedom = 4N 2 — (N + 
l) 2 — (N — l) 2 = 2A^ 2 — 2 which is found in a gauge theory on the fuzzy sphere without normal 
scalar field. 

The original U(N) gauge symmetry acts on the covariant derivatives D a as D 9 a = gD a g + , 
gE.U (N). This symmetry will be enlarged to the following U(2N) symmetry. First we introduce 
another covariant derivative D a = L a + A' a through a gauged idempotent 7 D ' given by a similar 
equation to (|32|). As before we will also set $ = 0. From the two idempotents jd and j d i we 
construct the link variable W as follows 

w = lD n D . (33) 

The extended U(2N) symmetry will then act on W as follows W — >VWV + , VeU(2N). This 
transformation property of W can only be obtained if we impose the following transformation 
properties j D ' — >V r y D >V + and 7d — >V^dV + on j D > and 7/5 respectively. Hence the U(N) 
subgroup of this U(2N) symmetry which will act on D a as D a — >gD a g + will also have to act 
on D a as D' a — >gD' a g + . Under these transformations the gauge fields A a and A' a transform as 
A a — >gA a g + + g[L a , g + ] and A' a — >gX a g + + g[L a , g + ] respectively like we want. Remark also 
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that for every fixed configuration A' a the link variable W contains the same degrees of freedom 
contained in j D . 

The main idea is that we want to reparametrize the gauge field on S 2 N in terms of the 
fuzzy link variable W and the normal scalar field $. In other words we want to replace the 
triplet (Ai, A 2 , A 3 ) with (W, $) where W is the object which contains the degrees of freedom 
of the gauge field which are tangent to the sphere. Thus in summary we have the coordinate 
transformation (A\, A2, A3) — >{W, $). We can check that we have the correct measure, viz 

J dA 1 dA 2 dA 3 oc J dWd$. (34) 

It remains now to show that the enlarged U(2N) symmetry reduces to its U (N) subgroup in the 
large N limit. The starting point is the 2N— dimensional one-plaquette actions with positive 
coupling constants A and A', viz 

S P = ^Tr 2N (W + W + -2) , S' P = -^Tr 2N (W 2 + W +2 -2). (35) 
A A 

We have the path integral 



Z P oc f rf 7D /rf$'5($') f dWd$5($)e s r +s 'r. (36) 

The extra integrations over 7^ and $' ( in other words over D a ) is included in order to 
maintain gauge invariance of the path integral. The integration over W is done along the orbit 
W = 7£>'7d inside the full U(2N) gauge group. In the large N limit this path integral can be 
written as 

Z P = f dA a 5(<$>') I dA a 5(<$>)e Sp+s 'r . (37) 

We need now to check what happens to the actions S P and S' P in the large iV limit. We 
introduce the 6 matrices 2A a = A a — A' a and 2A a = A a + A' a with the transformation laws 
A a — >gA a g + and A a — >gA a g + + g[L a ,g + ]. For the continuum limit of the action S P + S' P we 
obtain after a long calculation the effective theory 2 [T2] 

Z' P = J dAaSfax^ASje 3 ? (38) 

where 

Sf = AT 2 log(^)-^Tr(*[L a ^ 

1 (39) 

2 The path integral over the three matrices A a is dominated in the large N limit by the configurations A a = 0. 
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The coupling constant Ai ( which is assumed positive in this classical theory for simplicity ) is 
defined in terms of A abd A' by 

114 

-5T5-F (40) 

Notice that this effective action is invariant not only under the trivial original gauge trans- 
formation law A a — >A a but also it is invariant under the non-trivial gauge transformation 
A a — >A a + g[L a , g + ] where g&U(N). This emergent new gauge transformation of A a is iden- 
tical to the transformation property of a U(l) gauge field on the sphere. Therefore the action 
Sp given by the above equation is essentially the same U(l) action — (S — So) obtained from 
(JJJ) provided we make the following identification 

16 16 , 1 4 x 1 a A a 4 

(41) 



iV 2 Ai N 2 A A Ag 2 4 4iV 2 

between the U(l) gauge coupling constant g on the fuzzy sphere and the one-plaquette model 
coupling constant Ai- Let us also remark that in this large N limit in which g is kept fixed 
the one-plaquette coupling constant Ai goes to zero. Hence the fuzzy sphere action with fixed 
coupling constant g corresponds in this particular limit to the one-plaquette gauge field in the 
weak regime and agreement between the two is expected only for weak couplings ( large values 
of a ). 



5.2 The one-plaquette path integral 

Let us decompose the 2Nx2N matrix W as follows 

In particular W\ = Wi is an (iV+1) X (iV + 1) matrix, W 2 = W 2 is an (N — 1) x (iV — 1) matrix 
and W\2 is an (iV + 1) x (N— 1) matrix whereas the hermitian adjoint Wi 2 is an (N— 1) x (iV + 1) 
matrix. Since W + W = 1 we have the conditions 

W?Wi + W 12 Wf 2 = 1 , W 2 + W 2 + Wf 2 W 12 = 1 , WW12 + W 12 W 2 = 0. (43) 

Let us recall that since the integration over W is done along the orbit W = ^^^d inside U(2N) 
and since in the large iV limit both j D > and 7^ approach the usual chirality operator 7 = n a a a 
we see that W approaches the identity matrix in this limit. Thus we have the behaviour 
w i = {1d'1d)i — ►lw+i, ~W 2 = -(7d'7d)2 — ►Ijv-i and W 12 = (j D 'j D )i2 — >0. 

The main approximation adopted in ^2] consisted in replacing the constraint W = Jd'Jd 
with the simpler constraint W — >1 2 n by taking the diagonal parts W± and —W 2 to be two 
arbitrary, i.e independent of j D ', unitary matrices which are very close to the identities 1n+i 
and ljv-i respectively while allowing the off-diagonal parts W12 and to go to zero. We 
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observe that by including only W\ and —W2 in this approximation we are including in the limit 
precisely the correct number of degrees of freedom tangent to the sphere, viz 2N 2 . Thus in this 
approximation the integrations over $, $' and 7 D ' decouple while the integrations over W12 
and W^2 are dominated by W12 — W12 = 0- There remains the two independent path integrals 
over Wi and —W2 which are clearly equal in the strict limit since the matrix dimension of W\ 
approaches the matrix dimension of —W2 for large N. Thus the path integral Z' P reduces to 



Z P oc [Z P (X,X)Y (44) 



where 



Z P {\, A') = J dWi exp |yTr(Wi + - 2) - yTr(W^ 2 + -2)}. 



(45) 



5.3 Saddle point solution 



The path integral of a 2— dimensional U(N) gauge theory in the axial gauge A\ = on a 
lattice with volume V and lattice spacing a is given by Z P (X, oo) v/a2 where Z P (X, 00) is the 
above partition function (|4*5Jl for A' = 00, i.e the partition function of the one-plaquette model 
Sp = ^Tr(W\ + Wi — 2). Formally the partition function Zp(X, X') v ^ a2 for any value of the 
coupling constant A' can be obtained by expanding the model Si + S[ around A = 00. Thus it 
is not difficult to observe that the one-plaquette action Sp + S' P does lead to a more complicated 
U (N) gauge theory in two dimensions. 

Therefore we can see that the partition function Z' p of a U(l) gauge field on the fuzzy 
sphere is proportional to the partition function of a generalized 2— dimensional U(N) gauge 
theory in the axial gauge A\ = on a lattice with two plaquettes. This doubling of plaquettes 
is reminiscent of the usual doubling of points in Connes standard model. We are therefore 
interested in the N— dimensional one-plaquette model 

Z P (X,X') = J dWexp[^-Tr(W + W + -2) - ^Tr(W 2 + W +2 -2)). (46) 

Let us recall that dW is the U(N) Haar measure. We can immediately diagonalize the link 
variable W by writing W = TDT + where T is some U(N) matrix and D is diagonal with 
elements equal to the eigenvalues exp(i8i) of W. In other words = 5ijexp(i9i). The 
integration over T can be done trivially and one ends up with the path integral 

Z P (X,X') = Jlit^i^- (47) 

The action S^ is given by 

Sn = ^Eeos^-|E^s2^ + ^Eln(sm^) 2 -f + f • (48) 

i i iytj 
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Since the link variable W tends to one in the large N — >oo limit we can conclude that all the 
angles tend to in this limit and thus we can consider instead of the full one-plaquette model 
action (pfSj) the small one-plaquette model action 



* = --^ + ^"^-^ + 0(0% (49) 



A 2 is given by 



2 2 1 

Xr-Xr + S- <50) 



For the consistency of the solution below the coupling constant Ai must be negative ( as opposed 
to the classical model where Ai was assumed positive ) and as a consequence the coupling 
constant A 2 is always positive. As it turns out most of the classical arguments of section 4.1 
will go through unchanged when Ai is taken negative. Thus in this present quantum theory of 
the model we will identify the effective one-plaquette action Sp with the fuzzy sphere action 
S — S ( which is to be compared with the classical identification — SjF = S — So) and hence 
we must make the following identification of the coupling constants 

16 1 a 4 

(51) 



iV 2 Ai 4g 2 AN 2 ' 

Furthermore it is quite obvious that the expansion will only be valid for small angles Q{ in 
the range — Let us also note that the action ()49|) can be obtained from the effective 

one-plaquette model 



Sp = T^Tr(W eS + W c 

2 ^ ^ 2N 



V*E«**f--Vff- (52) 

A 2 i A 2 

For small 9* s in the range -l<6f<l the total effective one-plaquette action becomes 

S$ = -^U0f) 2 + ^rE ^ 1 4 jj +0((6^n (53) 

2 i i^j 

The action (JSHj) must be identical to the action (}4T?j) and hence we must have 6^ ff = 28i and 
Af = 4A 2 . 

The saddle point solution of the action (|52j) must satisfy the equation of motion 

c\ -r neS aeS 

^mtf-feca-L^J-. (54) 

In the continuum large N limit we introduce a density of eigenvalues p{9) and the equation of 
motion becomes 
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-^j sin 9 eS = J dT cS p(r cS ) cot cff - Teff . (55) 

By using the expansion cot = 2 Y^=i ( sin n6 cos rar — cos nO sin ra~) we can solve this equa- 
tion quite easily in the strong-coupling phase ( large values of A2 ) and one finds the solution 

1 1 

p(°cs) = — + — TTif cos6> eff . (56) 



2tt ttA 



However it is obvious that this solution makes sense only where the density of eigenvalues is 
positive definite, i.e for A| ff such that 

^>0 &{\f)* = 2 <*\* 2 = 0.5. (57) 



2n 7rA 2 

In the continuum large N limit where a 4 is kept fixed instead of Ai we can see that j- scales 
with N 2 and as a consequence A2 = — Ai = j^z- Thus the critical value A2* = 0.5 leads to the 
critical value 

64 

at = — = 128 & a* = 3.36 (58) 
A2 

which is to be compared with the observed value 

a* = 3.35±0.25. (59) 

This strong-coupling solution (J56)) should certainly work for large enough values of A 2 . However 
this is not the regime we want. To find the solution for small values of A2 the only difference 
with the above analysis is that the range of the eigenvalues is now [—9*, instead of [— ir, +n] 
where 9* is an angle less than tc which is a function of A2. It is only in this regime of small 
A2 where the fuzzy sphere action with fixed coupling constant g is expected to correspond to 
the one-plaquette model. Indeed the fact that W — »1 in the large iV limit is equivalent to the 
statement that the one-plaquette model is in the very weak-coupling regime. In the strong- 
coupling region deviations become significant near the sphere-to-matrix transition point. 
In the 'very' weak-coupling regime the saddle point equation reduces to 

2#i 2 _ 1 

aT^S^ (60) 

This problem was easily solved using matrix theory techniques in [TJJj. See also [23]. In the 
large iV — >oo we find the density of eigenvalues 



p(9) = —p\ 2 -9*. (61) 
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It is obvious that this density of eigenvalues is only defined for angles 6 which are in the range 
— \ / 2\2<6<\ / 2\2- However the value of the critical angle 9* should be determined from the 
normalization condition J_* e d6p{9) = 1. This condition yields the value 



The fuzzy one-plaquette third order phase transition happens at the value of the coupling 
constant A2 = 0.5 where the eigenvalues e tdi fill half of the unit circles. This half is due to the 
fact that Of = 26 i. 

In [T2] we also computed the predictions coming from this model for the free energy and spe- 
cific heat. We found very good agreement between the fuzzy one-plaquette model and the data 
in the weak-coupling phase and even across the transition point to the strong-coupling phase 
until the matrix-to-sphere transition point where deviations become significant. In particular 
the specific heat is found to be equal to 1 in the fuzzy sphere-weak coupling phase of the gauge 
field which agrees with the observed value 1 seen in our Monte Carlo simulation. The value 1 
comes precisely because we have two plaquettes which approximate the noncommutative U(l) 
gauge field on the fuzzy sphere. 

6 Conclusion 

In this article we have determined to a large extent the phase diagram of noncommutative 
U(l) gauge theory in two dimensions using the fuzzy sphere as a non-perturbative regulator. 
The central tool we employed was Monte Carlo simulation and in particular the Metropolis 
algorithm. 

We have identified three distinct phases. 1) A matrix phase in the strong coupling regime. 
The order parameter < TrX a X a >= in this phase. 2) A fuzzy sphere phase at weak coupling 
with order parameter < TrX a X a > 7^0 and with constant specific heat. 3) A new strong 
coupling fuzzy sphere phase. Here the fluctuations are around a fuzzy sphere background, i.e 
< TrX a X a > 7^0, in addition the specific heat is non-constant in this phase. 

The transition between the weak and strong coupling fuzzy sphere phases is third order. 
The other two transitions appear to be first order. We have clear numerical evidence for a 
jump in the internal energy < S > between the matrix and weak coupling fuzzy sphere phase. 
The corresponding jump in < S > has become smaller or disappeared in the strong coupling 
fuzzy sphere to matrix phase transition. However the order parameter < TrX a X a > still jumps 
discontinuously. We observe that for the m = model the specific heat becomes constant in 
both the strong coupling matrix phase ( Cv = |iV 2 ) and the weak coupling fuzzy sphere phase 
( Cv = N 2 ). As the mass m? increases a new third phase opens up and the three phases meet 
at a triple point. 

In this article we have also confirmed the theoretical one-loop prediction of the S 2 N — to- 
matrix critical line [22] • The transition between strong and weak couplings fuzzy sphere phases 
is found to agree with the -k expansion prediction of the one-plaquette critical line in the infinite 




(62) 
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mass limit. It seems that near these lines these approximations ( the one-loop and the one- 
plaquette ) are essentially exact. We also gave a Monte Carlo measurement of the triple point 
where the three phases meet. 

We would like to indicate that a high precision measurement of the one-plaqutte critical line 
and the triple point would be highly desirable. We also lack a theoretical control of the triple 
point. Improvement of the one-plaquette approximation of the NC U(l) gauge field on the 
fuzzy sphere S 2 N is necessary. In particular it would be very interesting to have an alternative 
more rigorous derivation of the one-plaquette critical value 3.35. Furthermore we believe that 
an extension of this approximation to 1) the regime of small masses and 2) the strong-coupling 
phase of the gauge theory is possible and needed. The Monte Carlo measurement and the 
one-loop theoretical description of the S%— to- matrix critical line are on the other hand very 
satisfactory. 

The most natural generalization of this work is Monte Carlo simulation of fuzzy fermions 
in two dimensions 27] and fuzzy topological excitations [26J. In particular our current project 
consists of the simulation of the NC Schwinger model and the NC two dimensional QCD on 
the fuzzy sphere. Then one must contemplate going to NC 4 dimensions with full QCD. Early 
steps towards these goals were taken in j2H| arid in the first reference of [15] . Supersymmetric 
models are also possible and in some sense natural |2H|. It would be nice to have Monte Carlo 
control over such supersymmetry. 
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A The order parameters and probability distribution 

In most of this section we restrict our discussion to the case m — 0. The case m^O has 
the same order parameters and probability distribution and we can show that they behave in 
exactly the same way. 

A.l Order parameters 

The model (J7J) is symmetric under U (N) gauge transformations of the matrices X a and as 
a consequence we can only attach a physical meaning to gauge invariant quantities which are 
constructed out of X a . In other words we have to measure gauge invariant observables. Let us 
introduce the scalar field $ defined by 



3 
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This field $ can be decomposed in the basis of NxN polarization tensors Y\ m as follows 



N-l I 

l=EE <i>lmY lm . (64) 
1=0 m=-l 

We remark that Y 00 — 1jv> ^i±i — ^k\/—L±, Yi = \/—L 3 and since X + = X we must also 

V _ y C2 y C2 

have (4>*)i m = (— l) m 0«-m- The total power in this field is given by 

I N-l I 

P= < v Tr¥ >=< £ E l^l 2 > • (65) 

iV 1=0 m=-l 

Another gauge invariant quantity we can measure is the power in the I = modes defined by 

P = < (^Tr$) 2 >=< fa > . (66) 
The data for P and Po is given in figure 12. The collapsed data is given in terms of P — 

C 2 

and Po = N2 i° as functions of a. From these results we can conclude that in the fuzzy sphere 
phase P = Po and thus the scalar field $ is proportional to the identity matrix since all its 
power is localized in the zero mode, i.e we have $ = ^^ =1 X n 2 = 0oo1jv- Furthermore a fit 
is given by P = P = a A c\ and hence we have essentially <poo — & 2 C2 in this phase which is 
consistent with the equilibrium configuration X a = aL a as expected. 
This result is confirmed by measuring the observables 

Pi =< jj TrX i > » etc - (67) 



The data for N = 6, 8 is shown on figure 13. The collapsed quantities are p\ = etc. We 

2 

find that in the fuzzy sphere phase we can fit the data to p a = ^y 2 - which is consistent with the 
number ±TrL\ = ±TrL\ = j^TrL\ = f. 

Let us now introduce the following 3 scalar fields ((a, b, c) = (1, 2, 3), (3, 1, 2), (2, 3, 1) ) 

^ ab =t[X a ,X b ] (68) 

The total powers associated with these scalar fields are given by 

j N-l I 

P ab = < -Tr$l b >=< £ £ |(0 afe U 2 > . (69) 

iV 1=0 m=-l 



In this case the powers in the I = modes vanish by construction. The definition of the modes 
' ' \b)im is obvious by analogy with equation ()64jl. The results are displayed on figure 14. The 

T 2 P C 
C2 

i til ;l ! I 1 o I I I ] (Hi Ti v> I if 

111 1 C OlW^li l v\~ <' V \ f > — 



collapsed quantities are P ab = — Pab . Again we can fit the data to the theoretical prediction 

4 

Pab = ^y 2 to a high degree of accuracy in the fuzzy sphere phase. Remark that the Yang- 
Mills action is given by < YM >= ^-{P\2 + P31 + -^23)- 111 the fuzzy sphere we clearly have 

<yM>=fp 12 = fp 3 i = f% 
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The powers Phat and Phat for N=6,8 



* 

* 
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Figure 12: The powers P, Pq for iV = 6, 8. 



We will also measure the following gauge invariant quantities 

Pcsi = -i< ^TrX 1 X 2 X 3 > , p cs2 = -i< ^TrX x X 3 X 2 > . (70) 
The results are shown on figure 15. In the fuzzy sphere phase we expect that the power p csl 

3 3 

behaves as p cs \ = whereas the power p cs2 behaves as p cs2 = — 9L § L - These are precisely 
the correct fits in the fuzzy sphere phase found for N = 4, 6 and 8 respectively. The collapsed 
powers are p cs i = ^^^ csl , etc. We remark that the Chern-Simons-like action is given by < 
CS >= —2aN 2 (p cs i —p cs2 ). In the fuzzy sphere phase we clearly have < CS >= —4aN 2 p cs i = 
AaN 2 p cs2 . 



A. 2 Geometric interpretation 

The covariant derivatives X a can in general be expanded in terms of NxN spherical har- 
monics Y\ m as follows 

3 oo I 

X a = X a L b + X a , X a = {x a )oQ%Q + S { X a)lmYl m - (71) 

6=1 1=2 m=-l 

The three vectors x a are the modes of X a with angular momentum I = 1 since Y%±i = 
±tJ~^L± = ±^J^(Li±.iL 2 ) and Y w = ^J^L 3 . They are vectors in R 3 which define the 
geometry of a parallelepiped. This geometry is precisely determined by the dynamics of X a 
given by the action (jZJ. The order parameters p a , P a b and p cs \ — p cs2 have the simple interpre- 
tation of the lenghts squared, the areas of the faces squared and the volume respectively of this 
parallelepiped. We have ( by setting X a = ) the following expressions 

i = pi = -a x l , etc 
3 
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The powers phat a for N=6,8 
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Figure 13: The powers p a for N = 6,8. 



a 2 = Pi 2 = -a 4 (xixi 2 ) 2 , etc 
3 

v =pcsi ~Pcs2 = ^a 3 fi.(f 2 xf 3 )- 

(72) 

The full effective action in terms of the vectors x a takes the form 

~4 

c 2 a 



S[x a ] 



- 4r l_ 1_ 1 



-(fiXi 2 ) 2 + -(iiXi 3 ) 2 + -(f 2 xi 3 ) 2 — 2xi.(x 2 y<x 3 ) — m 2 x\ — m 2 x\ — m 2 x\ 



3 

— r~ ~ ^a/jcd(^i^i ~~ t"~ *^'2*^'2 ~~ ^" ^3*^3) v'^'l^'l ~~ ^" *^2*^2 ~~ ^" ^3*^3) 



2 

m 



+ S[x a }. (73) 



5[^ a ] is the quantum action obtained by integrating out the field X a ( equation (fTTj) ) from 
the theory. The coefficients d a b c d can be computed easily from the definition j^TrL a LbL c Ld = 

2 

^■dabcd- The original U(N) gauge symmetries are now implemented by 0(3) orthogonal symme- 
tries which take the 3— dimensional vectors x a to = Rx a . The effect of quantum fluctuations 
in this problem is to deform the shape of the parallelepiped. In particular the first order phase 
transition from the fuzzy sphere to the matrix phase is now seen as the transition where the 
parallelepiped collapses. In terms of the lengths squared I 2 , the areas of the faces squared a 2 
and the volume v we have 



I 2 = 1 < TrX 2 >= I | fuZZy Sphe ^ Phase \ , etc. (74) 
c 2 [Im matrix phase. 



2 ^ , T [y yl 2 v y fuzzy sphere phase i 

a = <Tr[Xi,X 2 \ >= , . , } ,etc. (75) 

C2 I a m matrix phase. 1 



33 



The powers Phat ab for N=6 and £ 
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Figure 14: The powers P a b for N = 6,8. 



v = < Tr Xl [X 2) X 3 ] >= { * fUZZy Sphe ^ Phase ) . (76) 

4c 2 { v m matrix phase. J 

From the data we can see that the areas of the faces squared o? m and the volume v m in the 
matrix phase are constant approaching the values 2 and respectively for small values of the 
coupling constant a. However the length squared l 2 m scales as N~? and thus it becomes in 
the limit. 



A. 3 Probability distribution 

As we said before for m = we can take X a to be traceless without any loss of generality 
and consider only the probability distribution and the partition function given by 

V[X a ] = 6(TrX ^ S[Xa] , m = J [dX a ]5(TrX a )e- s ^. (77) 

The classical absolute minimum of the model is given by X a = aL a . The quantum minimum is 
given by X a = acpL a where acfi plays the role of the radius of the sphere with a classical value 
equal a. The complete one-loop effective potential in this configuration is given in the large N 
limit by the formula (jlOj) . The solution (f> of the equation of motion 4 — (f> 3 + -3$ = approaches 
the classsical value 1 as one increases the coupling constant a much above the critical value a*. 
Indeed it is not difficult to check that up to the order of 4% we have 

2 12 1 

<P = l---- + 0(—). (78) 
a 4 a 8 a Vi 
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In this section we report on the measurement of the radius of the fuzzy sphere. A natural 
definition of the radius of is given by the observable 

* = ^<^TrXl>. (79) 

The aim now is to make a precise measurement of by measuring R 2 and its probability 
distribution V(R 2 ). Numerically we thermalize and then we take Tmont measurements of R 2 , 
we determine the minimum and maximum values R 2 mi and R 2 ma respectively and divide the 
interval [R^, R^a] into q = 2 6 + 1 smaller intervals of equal length 5 = Rma q Rmi . For every 
measurement R? i = 1, Tmont, we compute the integer 



3 



( R; — R 



integer part I 



'■nil 



It is clear that the value R 2 will lie exactly in the j—th interval, in other words 

R 2 = R 2 mt +jS. (81) 

We count the number of times N(J) we get the value R\ and we define the corresponding 
probability V(j) by 

V(j) = 7^^. (82) 

KJJ Tmont v ; 

Remark that this probability satisfies Y?j=o'P(J) = 1- -'■ n other words for all j = 0, q we have 

r(j)<i. 

We observe two radically different behaviour depending on wether we are inside the fuzzy 
sphere phase or the matrix phase. Figure 16 shows the probability distribution in the fuzzy 
sphere phase whereas figure 17 shows the probability distribution in the matrix phase. 
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Once again we find a good agreement between the theory and the simulation in the fuzzy 
sphere phase. More precisely we find that the value of R 2 depends explicitly but slowly on the 
coupling constant a but it does not depend on N. In figure 13 we see clearly that the value 
of R 2 at the peak of the probability V is increasing with increasing a. We also observe clearly 
how the value of R 2 at the peak is more or less the same for a given value of at with different 
N. These results are consistent with equation (|78jl. 

In figure 17 we plot the probability distribution V as a function of r 2 = a 2 C2R 2 for N = 6, 8 
and for values of the coupling constant which are less than 2. In other words we are inside the 
matrix phase. For a = 0.5 , 1 and 1.5 and for all the values of N we observe that the probability 
distribution in this phase peaks essentially around the same value which is estimated to be in 
the range r 2 = 2.4 — 2.8. Hence we can conclude immediately that for a fixed value of the 
coupling constant a inside the matrix phase the order parameter R 2 = ^4^- will be peaked 
around smaller and smaller values as we increase N. This means in particular that the Chern- 
Simons-like term in the action is playing no role in this matrix phase and as a consequence we 
have no an underlying spacetime structure of a fuzzy sphere. 

The results are summarized as follows: 

i? 2 = i<iV TrX 2 >= { 1 fuZZy Sphe ^ Phase ) . (83) 
a 2 c 2 N ^ [0 matrix phase. J v ' 
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